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ABSTRACT 

We present new determinations of effective temperature, surface gravity, and masses for a sample of 
46 hot DA white dwarfs selected from the EUVE and ROSAT Wide Field Camera bright source lists in 
the course of a near-IR survey for low mass companions. Our analysis, based on hydrogen NLTE model 
atmospheres, provides a map of LTE correction vectors, which allow a thorough comparison with previous 
LTE studies. We find previous studies underestimate both the systematic errors and the observational 
scatter in the determination of white dwarf parameters obtained via fits to model atmospheres. 

The structure of very hot or low mass white dwarfs depends sensitively on their history. To compute 
white dwarf masses, we thus use theoretical mass-radius relations that take into account the complete 
evolution from the main sequence. We find a peak mass of our white dwarf sample of O.SQAf©, in 
agreement with the results of previous analyses. However, we do not confirm a trend of peak mass with 
temperature reported in two previous analyses. 

Analogous to other EUV selected samples, we note a lack of low mass white dwarfs, and a large 
fraction of massive white dwarfs. Only one white dwarf is likely to have a helium core. While the lack of 
helium white dwarfs in our sample can be easily understood from their high cooling rate and therefore 
low detection probability in our temperature range, this is not enough to explain the large fraction of 
massive white dwarfs. This feature very likely results from a decreased relative sample volume for low 
mass white dwarfs caused by interstellar absorption in EUV selected samples. 

Subject headings: stars: atmospheres — stars: evolution — stars: white dwarfs — binaries: close 



1. INTRODUCTION 

UV observations of EUV-detected stars have revealed 
the presence of about 15 hot white dwarf companions 
to bright stars in non-interacting binary systems (e.g., 
Burleigh, Barstow, & Fleming 1997). At optical wave- 
lengths, these WDs arc hidden because of their close prox- 
imity to much more luminous companions which are main 
sequence (spectral type K or earlier) or evolved stars. 

A fascinating variety of objects are known or proposed 
to contain white dwarf (WD) stars in interacting binary 
systems. A partial list includes novae, cataclysmic vari- 
ables, symbiotic stars, Ba and CH giants, Feige 24-type 
systems and dwarf carbon stars (Green & Margon 1994). 
These systems offer great insights to evolution and dynam- 
ical processes in binaries. 

A number of interacting binary systems where the 
white dwarf is the primary (i.e., optically brightest) star 
have also been found among EUV-detected systems (e.g., 
6 close, interacting white dwarf/red dwarf binaries by 
Venues & Thorstensen 1994). Optical or ultraviolet spec- 



tral observations are most commonly used to detect com- 
panions to WD primaries, by searching for (1) the presence 
of narrow Balmer line emission overlying the broad smooth 
Balmer absorption of the WD, (2) a composite white dwarf 
-|- main sequence spectrum, or (3) radial velocity (RV) 
variations. However, only WDs with very close, or intrin- 
sically active companions will be found by method ( 1 ) . For 
hot white dwarf systems, composite spectra (2) are only 
expected to be visible if the companion's spectral type is 
early enough. RV variations (3) require multiple observa- 
tions at high spectral resolution, and detection strongly 
favors close and/or massive companions. 

All the discoveries mentioned above have been strongly 
dominated by these selection effects, with companions bi- 
ased to earlier types than predicted by the simulations of 
deKool & Ritter (1993) and others. Scaling from deKool 
& Ritter's results, Vennes & Thorstensen (1994) estimate 
that "at least twice as many close binary systems remain to 
be identified from EUV surveys, most of them with a low 
mass secondary." The resulting sample of binaries known 
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to date therefore must diverge strongly from the intrinsic 
distribution, in overall normalization, as well as in mass 
and spectral type of the main sequence companions. 

The current study, conceived as a complement to op- 
tical studies, began as a near-IR photometric survey for 
low mass companions to hot white dwarfs (WDs). By 
investigating only EUV-detected WDs, we obtain a very 
reasonably-sized but complete sample of young WDs, next 
to which very late-typo dwarf companions can be detected 
in the near-infrared by searching for a K excess over that 
expected from the white dwarf. Many hot white dwarfs 
(Teff > 24000K; Finley et al. 1993) have been detected in 
the recent EUV all-sky surveys. EUV detection of these 
hot WDs depends primarily on their temperature, dis- 
tance, and the intervening Galactic ISM. Our sample of 
EUV WDs (whose selection we define below) offers excel- 
lent flux contrast in the IR relative to optical; cool com- 
panions will almost always be brighter in the K band than 
the hot WDs. 

To know what K mag to expect for the WDs, we benefit 
from constraints on log 3, radius, and T^s derivable from 
optical spectra for the WDs in our sample, using NLTE 
model atmosphere fits (Napiwotzki et al. 1993, Napiwotzki 
1997). The resulting predictions for K magnitudes allow a 
direct search for any IR excess from a cool companion. In 
some cases, IR colors will also provide a preliminary spec- 
tral type. Results from the IR survey will be presented in 
an upcoming paper. 

An additional motivation is the study of the white dwarf 
mass distribution. Since the pioneering work of Koester, 
Schulz & Weidemann (1979; KSW) it is well estabhshed 
that the masses of white dwarfs cluster in a narrow range 
around Q.QMq, remarkable given that white dwarfs stem 
from progenitors with masses ranging from below IMq 
up to f«8M0. Precise knowledge of the white dwarf mass 
distribution puts constraints on the theory of stellar evo- 
lution, especially the poorly understood mass loss process 
during the final stages of stellar evolution. With two recent 
exceptions (Beauchamp et al. 1996 and Dreizler & Werner 
1996, who analyzed samples of helium-rich DB and DO 
white dwarfs, respectively) the mass distribution has only 
been determined for hydrogen-rich DA white dwarfs. How- 
ever, this is not a severe limitation, because this spectral 
class comprises about 80 % of all known white dwarfs. 

The analysis of KSW, along with other follow-up in- 
vestigations in the early eighties, used photometric data 
of which the Stromgren and Greenstein multichannel 
colors were the most suitable. Both systems provide 
temperature- and gravity-sensitive indices. Alternatively, 
KSW and others used trigonometric parallax measure- 
ments to directly calculate the stellar radius. However, the 
latter method is practical only for a small sample and suf- 
fers from considerable measurement uncertainty. Unfor- 
tunately, the photometric indices have their highest sensi- 
tivity near 10000 K. At this temperature DA white dwarfs 
have a convectivc atmosphere and the resiilts depend crit- 
ically on the adopted parameters of mixing length theory. 

The situation improved at the beginning of the nineties, 
when the development of modern, highly efficient detec- 
tors made it possible to obtain high-quality spectra of large 
numbers of white dwarfs and determine the stellar param- 
eters from a fit to the detailed profiles of the Balmer lines. 



This method yields sufficient accuracy for white dwarfs hot 
cnoiigh to have a radiative envelope. The first comprehen- 
sive sample of white dwarfs analyzed by this method was 
presented by Bergeron, Saffer & Liebert (1992; hereafter 
BSL). As attributed to higher precision of spectroscopic 
methods, this investigation yielded a white dwarf mass 
distribution even narrower than found by KSW and other 
previous studies. 

At the same time, Kidder (1991) analyzed a sample of 
hot DA white dwarfs discovered through positional coinci- 
dences of catalogued hot DA white dwarfs in existing soft 
X-ray databases. Three soft x-ray sources corresponding 
to white dwarfs were found, having relatively low effec- 
tive temperatures, f»25,000K, determined independently 
using complementary optical and UV spectroscopy. Kid- 
der et al. (1992) analyzed an expanded sample to derive 
photospheric He abundances for the hotter objects and es- 
tablish an effective observational low-temperature thresh- 
old for the detection of pure hydrogen DA white dwarfs at 
soft X-ray wavelengths. 

In 1997 three groups (Marsh et al. 1997, M97; Vennes 
et al. 1997, V97; Finley, Koester & Basri 1997, FKB) pub- 
lished results on the mass distribution of Extreme Ultra- 
violet (EUV) selected white dwarfs. Due to the selection 
criterion, these samples contain the hottest white dwarfs 
(Teff> 25000 K), as cooler white dwarfs do not emit sig- 
nificant EUV radiation. The derived mass distributions 
in the EUV-selectcd samples are similar to that of BSL, 
but show some interesting deviations in detail. The fre- 
quency of very high mass white dwarfs is much larger, 
and that of very low mass white dwarfs much smaller, 
than in BSL. These findings can at least partly be ex- 
plained by selection effects (sec the discussion in FKB). 
More serious is a trend of the peak mass with temperature. 
V97 found that their mass distribution peaks at O.598M0, 
while the BSL distribution peaks at O.568M0, with masses 
computed using Wood's (1995) mass-radius relation with 
"thick" layers (Afn = IO^^'Mwd, Afnc = lO^^Afwo)- 
This discrepancy diminishes slightly, if the "very thin 
layer" (Mne = 10~^Mwd, no hydrogen layer) mass-radius 
relations are used (peak masses of O.556M0 and O.532M0 
for the V97 and BSL sample, respectively). V97 inter- 
preted this as evidcncx; for a very thin hydrogen layer of 
the DA white dwarfs. However, the effects are small so 
this result depends strongly on the accuracy of the derived 
stellar parameters. 

FKB estimated the internal accuracy of different anal- 
ysis methods from Monte Carlo simulations. The preci- 
sion reachable by Balmer line fitting is very compelling: 
Arcff/Tcff< 0.01 for Teff< 60000 K. However, for spectra 
with very high signal-to- noise ratios (S/N), errors intro- 
duced by details of the observation and reduction tech- 
niques (e.g., extraction, flat flelding, flux and wavelength 
calibration) might be more important, but are very diffi- 
cult to determine. Additionally, one has to take into ac- 
(;ount differences in the model atmosphere caknilatioiis and 
fitting procedure. Together with the results presented in 
this paper, we now have four samples of hot white dwarfs, 
analyzed in a similar way, and with significant overlap. 
This offers the opportunity to determine the real accu- 
racy of the spectral analysis of hot white dwarfs, including 
many possible systematic effects. 



3 



We present the selection criteria of our sample in Sect. 2 
and the observations and data reduction procedures in 
Sect. 3. Details on our model atmospheres are given in 
Sect. 4. The results and a detailed comparison with the 
previous analyses of EUV selected white dwarfs are pre- 
sented in Sect. 5. We finish with a discussion of our results 
and an outlook. 



2. SAMPLE SELECTION CRITERIA 

We chose to limit our uniform sample to DAs, for which 
models provide the best temperature and mass constraints. 
We start with 73 known DA white dwarfs in the EUVE 
bright source list (Malina et al. 1994). Excluding sources 
at low galactic latitudes < 15) and southerly dechna- 
tions {S < —20) yields a list of 28 DAs. A similar proce- 
dure for non-overlapping DAs listed in the ROSAT Wide 
Field Camera survey Bright Source Catalogue (Pounds ct 
al. 1993) yields 29 objects. 

We have removed from our uniform sample 2 well- 
known stars with published sensitive optical spectropho- 
tometry and IR photometry (Fcigc 24 and HZ 43). Two 
DAs with broad line profiles due to magnetic splitting 
were also excluded (PC 1658-1-441 and PG0136-h251). 
Eight known binaries are also excluded from the uniform 
sample: V471 Tau (Vennes, Christian, & Thorstensen 
1998), PC 0824+289, (Heber et al. 1993) HD74389B 
(Liebert, Bergeron, & Saffer 1990), REJ1016-052 (V97), 
PC 1033-^464 (CD 123), (Green, Schmidt, & Liebert 
1986), REJ1426-h500 (V97), REJ1629-h780 (Catalan et 
al. 1995), and IK Peg (Wonnacott et al. 1993), leaving 47 
objects. 

In this paper, a handful of objects which fell out- 
side the uniform sample definitions just outlined were in- 
cluded for observation. These include the known binaries 
PC 1033+464 and RE J1629+780 and the magnetic white 
dwarf PC 1658+441, as well as MCT 0455-2812 which was 
outside the sample declination limits. 

Due to observing constraints (a combination of weather, 
poor seeing and faint objects, or celestial placement 
of objects) no spectra were obtained for sample ob- 
jects REJ0443-034, REJ0916-194, or PG1040+451. 
PG1234+482 was originally classified as a sdB star and 
thus excluded; it has since been reclassified as a DA (Jor- 
dan et al. 1991). The final sample we analyze here thus 
includes 46 DA white dwarf stars for which we present new 
model NLTE fits to optical spectra. 

We note that since the sample selection was performed, 
several relevant discoveries pertaining to sample objects 
have been made. REJ0134-160 (GD984) has central 
Balmer emission components produced by a dMe compan- 
ion. RE J1440+750 turns out to be a magnetic DA (Drei- 
zler et al. 1994). The uniform sample will be discussed in 
a foUowup paper treating the IR photometry and binary 
fraction. 



3. OBSERVATIONS 

Dates are listed for all observations in Table |l|. On the 
nights of 04 - 06 January 1996 we obtained spectra at Stew- 
ard Observatory's Kitt Peak Station using the Bok 2.3-m 
reflector equipped with the Boiler & Chivens Cassegrain 
spectrograph and UV-flooded Loral 800x1200 CCD. Most 



spectra were dispersed with a 600 I /mm first-order grat- 
ing used behind a 4.5" x 4' long slit. The instrumentation 
provided wavelength coverage AA3400-5600 at a spectral 
resolution of ^5A FWHM. On the last night of the observ- 
ing run we employed a new 400 1/mm grating providing 
wavelength coverage AA3500-6790 at a spectral resolution 
of -71 FWHM. 

We also obtained spectra at the MMT on Mt Hopkins, 
09 - 11 May 1996, using the 300 I /mm grating in first 
order on the Loral 3kxlk CCD of the Blue Channel spec- 
trograph. This yields coverage from about 3350 to 8800 A, 
and the 2" slit width we used resulted in a spectral reso- 
lution of about 4 A FWHM. Several objects were kindly 
obtained for us using the identical instrumental configura- 
tion at the MMT by Perry Berlind on 08 April 1997. 

Exposures at both telescopes ranged from one to thirty 
minutes for program stars, and for all observations the 
long slit was rotated to the parallactic angle according 
to the calculations of Filippenko (1982). The airmasses 
were held below 1.5 in almost all cases. All spectra were 
extracted from the two-dimensional images and reduced 
to linear wavelength and intensity scales using standard 
reduction packages in the Image Reduction and Analysis 
Facility (IRAF). These operations included bias subtrac- 
tion, flat field division by images obtained by exposing on 
dome or internal quartz lamps, centroiding and summa- 
tion of the stellar traces on the two-dimensional images, 
sky subtraction, wavelength calibration using spectra of 
He/Ar arc lamps, and absolute flux calibration using the 
flux standards of Massey et al. (1988). Details of these 
reduction procedures are given by BSL. 

We present in Figs. |l and ||our collection of white dwarf 
spectra. The spectra obtained at the Steward observatory 
2.3 m telescope are found in Fig. |^, and the MMT spectra 
in Fig. Ij. Some stars observed repeatedly appear in both 
figures. The signal-to-noise ratio of our spectra ranges 
from 35 to 200 with an average of 90. Ninety percent of 
our spectra have a S/N of at least 45. 

4. MODEL ATMOSPHERES 

We calculate hydrogen model atmospheres with the 
NLTE code developed by Werner (1986). Basic assump- 
tions are those of static, plane-parallel atmospheres in hy- 
drostatic and radiative equilibrium. In contrast to the at- 
mospheres commonly used to analyze DA white dwarfs, we 
relax the assumption of local thermal equilibrium (LTE) 
and solve the detailed statistical equilibrium instead. As 
described in Werner (1986), the accelerated lambda iter- 
ation (ALI) method is used to solve the set of non-linear 
equations. The impact of NLTE on white dwarf atmo- 
spheres is discussed in detail in Napiwotzki (1997). 

H I levels and lines are included in NLTE up to n = 16. 
Line blanketing by the Stark broadened hydrogen lines is 
taken into account consistently. As the hydrogen atmo- 
spheres of DA white dwarfs are stable for Tcg> 15000 K, 
convection is not included in our atmospheric models. 
Pressure dissolution of the higher levels is described by the 
Hummer & Mihalas (1988) occupation probability formal- 
ism following the NLTE implementation by Hubeny et al. 
(1994). The synthetic spectra are computed with the ex- 
tended VCS broadening tables (Vidal et al. 1970) provided 
by Schoning & Butler (priv. comm.) and Lemke (1997). 
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We followed the prescription of Bergeron (1993) and in- 
creased the critical ionizing field adopted to calculate the 
occupation probability by a factor two. The motivation 
is not a flaw in the Hummer & Mihalas (1988) formalism, 
but a compensation for the inadequacy of the standard 
Stark broadening theory when line wings overlap. Our 
NLTE model grid covers the temperature range 17000 K < 
Teff<10G0G0K (stepsize increasing with Tcff from 2000 K 
to 10000 K) and gravity range 6.50 < logg < 9.75 (stepsize 
0.25). 

Although deviations from LTE are small for most DA 
white dwarfs they become significant for the hottest stars 
in our sample (cf. Napiwotzki 1997). Since we intend 
to compare our results with three other samples analyzed 
by means of LTE atmospheres, we have produced a map 
with LTE correction vectors. For this purpose we calcu- 
lated a set of LTE model atmospheres using the technique 
described in Napiwotzki (1997) of drastically enhancing 
collisional rates between the atomic levels in the NLTE 
code. This forces the occupation numbers to be in LTE 
and guarantees consistency with the NLTE atmospheres. 

The synthetic LTE spectra were transformed into "ob- 
served" spectra by convolving them with a Gaussian of 5 A 
FWHM, rebinning them to 2 A and adding Poisson noise 
corresponding to a continuum S/N of 100. These simu- 
lated spectra were analyzed with the NLTE grid following 
the procedure outlined in Sect. |^. 1000 simulations were 
run for every parameter set to eliminate the effect of ran- 
dom errors. The resulting offsets are displayed in Fig. |[ 
The orientation of the vectors corresponds to the correc- 
tion, which must be applied to transform LTE results to 
the NLTE scale. 

As expected from the results of Napiwotzki (1997) the 
differences are negligible for DA white dwarfs cooler than 
«40000K, but can be significant for hotter stars. The 
LTE vectors show the trend that NLTE effects increase 
with increasing temperature and decreasing gravity (see 
e.g. Fig. 3 and 4 in Napiwotzki, 1997, for the case of DA 
white dwarfs). The small corrections found for the mod- 
els with the highest temperature and lowest gravity seems 
to contradict this behavior. A look at the line profiles 
reveals that the NLTE deviations are larger for, say, the 
Toff= 90000 K, log 5 = 7.0 model than the 70000 K, 8.0 
model. However, the correction vectors are not a simple 
function of NLTE deviations measured e.g. as equivalent 
width difference, but depend also on the way line profiles 
vary with temperature and gravity. The small corrections 
found for the high temperature/low gravity models are 
produced by the cancellation of these effects. We expect a 
reliable transformation between the LTE and NLTE tem- 
perature scales only if the deviations are not too large. 
A conservative upper limit is R:!70000K. A thorough com- 
parison with the previous LTE analyses is presented in 
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We checked our models by a comparison of our LTE 
spectra with some DA model spectra kindly provided by 
D. Koester. Model parameters were T^s = 30000 K to 
70000 K in 10000 K steps and logg = 8.0 for pure hydro- 
gen models. Input physics are very similar. In particu- 
lar both model calculations adopt twice the critical ion- 
izing field for the calculation of line profiles. We treated 
Koester's model spectra the same way we treated the LTE 



spectra above and fitted them with our LTE grid. The 
result was quite satisfactory: the temperature differences 
were always below 1.5% and the gravity differences never 
exceeded 0.03 dex. 

Metals were ignored in our calculations, but they can 
modify the hydrogen line profiles by their effect on the 
atmospheric structure. Lanz et al. (1996) analyzed the 
Balmer lines of the hot DA G 191 B2B with pure hydrogen 
LTE and NLTE atmospheres and an NLTE model with full 
metal line blanketing. They concluded that the effect of 
metal line blanketing on the Balmer lines was relatively 
small, and the difference between LTE and NLTE was 
found to be the most important effect. A recent study 
by Barstow et al. (1998), which investigated several hot 
white dwarfs in the temperature range around 60000 K 
derived larger metal line blanketing effects of the order of 
the NLTE effects. Since the LTE analyses of M97, V97, 
and FKB are based on pure hydrogen models, our results 
should be consistent with theirs in any case. 

5. SPECTRAL ANALYSIS AND RESULTS 

Atmospheric parameters of our DA white dwarfs are ob- 
tained by simultaneously fitting line profiles of the ob- 
served Balmer lines with the NLTE model spectra de- 
scribed above. We use the least-square algorithm de- 
scribed in BSL. The observed and theoretical Balmer line 
profiles are normalized to a linear continuum (both spec- 
tra are Fx) in a consistent manner. Wavelength shifts are 
determined with a cross-correlation method and applied 
consistently to each complete spectrum. The synthetic 
spectra are convolved to the observational resolution with 
a Gaussian and interpolated to the actual parameters with 
bicubic splines, and interpolated to the observed wave- 
length scale. 

The atmospheric parameters Tcs and logg are then 
determined by minimizing the value by means of a 
Levenberg-Marquardt steepest descent algorithm (Press et 
al. 1986). Several tests revealed that our interpolation rou- 
tine is rather robust concerning spacing of our model grid 
and yields reliable results even at the edge of the model 
grid. 

Finally, an estimate of the internal errors can be de- 
rived from the covariance matrix. In contrast to BSL, we 
estimate the noise of the spectra (ct) used for the fit 
from the neighboring continuum of each line. The S/N is 
adopted to be constant throughout the line. 
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Table 1 

ATMOSPHERIC PARAMETERS OF PROGRAM STARS 



tXtli J 


AT 

Other JNames 


-teflF 




l\/f / A/f 
IVl 1 IVIq 


Observation 
Date 


0007+331 


GD2 


46493+514 


7.83+0.05 


0.602 


6 Jan 1996 


0134-160'' 


GD 984, PHL 1043 


44866+667 


7.77+0.05 


0.572 


5 Jan 1996 


0237-122 


PHL 1400 


32077+177 


8.45+0.04 


0.890 


5 Jan 1996 


0348-005 


GD 50 


39508+464 


9.07+0.06 


1.215 


4 Jan 1996 


0427+740 




48587+1044 


7.93+0.08 


0.646 


5 Jan 1996 


0457-280 


MOT 0455-2812 


51199+786 


7.72+0.05 


0.570 


4 Jan 1996 


0512-004 




31733+139 


7.40+0.03 


0.458 


5 Jan 1996 


0521-102 




33186+301 


8.60+0.06 


0.980 


5 Jan 1996 


0841+032 




38293+252 


7.75+0.03 


0.544 


4 Jan 1996 


0902-040 




23218+160 


7.84+0.02 


0.544 


4 Jan 1996 


0907+505 


PG 0904+511 


32167+338 


8.11+0.07 


0.695 


4 Jan 1996 


0940+502 


PG 0937+506 


36034+283 


7.69+0.04 


0.519 


9 Apr 1997 


0957+852 




51311+1348 


8.37+0.10 


0.866 


5 Jan 1996 


1019-140 




31524+102 


7.92+0.02 


0.606 


9 Apr 1997 


1029+450 


PG 1026+454 


35518+247 


7.70+0.04 


0.520 


9 Apr 1997 


1032+532 




43587+506 


7.95+0.05 


0.644 


4 Jan 1996 


1033-114 


G 162-66, LTT3870 


24685+252 


7.85+0.03 


0.553 


4 Jan 1996 


1036+460'= 


GD123 


29361+251 


8.02+0.05 


0.648 


4 Jan 1996 


1043+490^= 




41132+1178 


7.94+0.13 


0.635 


5 Jan 1996 


1044+574 


PG 1041+580 


30338+153 


7.81+0.03 


0.550 


4 Jan 1996 


1100+713 


PG 1057+719 


41104+814 


7.84+0.09 


0.593 


4 Jan 1996 


1112+240 


Ton 61 


39824+636 


7.78+0.07 


0.563 


5 Jan 1996 


1122+434 


PG 1120+439 


26996+151 


8.31+0.02 


0.803 


9 Apr 1997 


1126+183^ 


PG 1123+189 


54334+1983 


7.76+0.13 


0.594 


4 Jan 1996 


1128-025 


PG 1125-026 


30699+380 


8.24+0.08 


0.767 


4 Jan 1996 


1148+183 


PG 1145+188 


25758+299 


7.91+0.04 


0.587 


4 Jan 1996 


1235+233^ 


PG 1232+238 


46569+523 


7.83+0.05 


0.602 


11+12 May 1996 


1257+220*^ 


GD153 


38926+142 


7.78+0.02 


0.560 


4 Jan 1996, 12 May 1996 


1336+694 


PG 1335+701 


29607+87 


8.34+0.02 


0.824 


10 May 1996 


1431+370 


GD 336 


34404+115 


7.91+0.02 


0.608 


11 May 1996 


1446+632 




37947+254 


7.79+0.04 


0.562 


11 May 1996 


1629+780^= 




41043+338 


7.92+0.04 


0.627 


11 May 1996 


1638+350 


PG 1636+351 


35404+142 


7.98+0.02 


0.642 


10 May 1996 


i04oH-41 i 




0001 rr 1 Qi 


OO-LH no 


U. / OO 


y May iyyo 


1650+403^ 




38144+211 


7.97+0.03 


0.643 


9 May 1996 


1711+664^ 




48989+757 


8.89+0.06 


1.141 


12 May 1996, 9 Apr 1997 


1726+583'^ 


PG 1725+586 


53561+542 


8.23+0.04 


0.795 


9 May 1996 


1800+683 


KUV 18004+6836 


44723+424 


7.80+0.04 


0.585 


10 May 1996 


1820+580 




44099+264 


7.78+0.03 


0.574 


9 May 1996 


1845+682 


KUV 18453+6819 


36120+189 


8.23+0.03 


0.770 


10 May 1996 


2116+735*^ 


KUV 21168+7338 


50812+354 


7.72+0.03 


0.569 


5 Jan 1996, 9 May 1996 


2207+252 




26964+174 


8.27+0.03 


0.779 


5 Jan 1996 


2312+104 


GD 246 


53088+968 


7.85+0.07 


0.624 


4 Jan 1996 



NOTES: The error estimates were taken from the procedure. 

'^spectrum sHghtly contaminated by cool companion; the cores of H/3 and H7 are excluded from the 

fit 

''Weighted mean of individual observations 

'^red part of spectrum contaminated by cool companion; H/3 and the core of H7 excluded from fit 

•^Close-by optical companion; spectrum apparently not contaminated 

*red part of spectrum contaminated by cool companion; H/3 excluded from fit 
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The results are given in Table ^, with illustrative exam- 
ples shown in Fig. ^ We adopt, for the moment, the usual 
practice and indicate in Table pj the internal errors esti- 
mated from the quality of the x fit. However, one should 
keep in mind that these errors can only serve as lower lim- 
its. We will show below (Sect. |5.3| ) that these formal errors 
derived from the fit significantly underestimate the real 
errors. 

External errors can be estimated from multiple obser- 
vations and analysis of the same star. We obtained repeat 
observations for a subsample of 6 stars, for which results 
are given in Table ||. The gravity values of all six stars 
agree within the estimated internal errors. The same is 
true for four temperature comparisons, too. However, the 
differences found for REJ1650-h403 and REJ2116-h735 
are significantly larger. This is in line with the external 
errors we estimate from a com pari son with the studies of 
M97, V97, and FKB (see Sect.lO). 



5.1. Binaries 

Five stars in our sample, RE J0134-160, RE J1036-I-460, 
REJ1043+490, REJ1126+183, and REJ1629+780, show 
clear signs of binarity in the red part of their spectra. Two 
more stars, REJ1711-h664 and RE J2207-I-502, are mem- 
bers of visual binaries. 

REJ1629+780: The red part of the spectrum of 
RE J1629-I-780 is heavily contaminated by a M type main 
sequence companion. The composite spectrum and the 
spectrum of the M star after subtracting the white dwarf 
component are shown in Fig. The characteristic bands 
of TiO are easily recognizable. Catalan et al. (1995) deter- 
mined a spectral type dM4. The Ha and H/3 lines are seen 
in emission, which indicates a chromospherically active Me 
star. Sion et al. (1995) detected a flare-like increase of the 
Balmer line emission. 

The comparison in Fig. |^ demonstrates that the blue 
part (A < 4400 A) of the white dwarf spectrum is not 
disturbed by the M dwarf. Thus, we excluded H/? and 
the core of H7 from the fit and derived the parame- 
ters reft= 41000 K and log 5 = 7.92 from H7 to H??. 
These parameters are in reasonable agreement with the 
results of Catalan et al. (Toff= 41800 K; \ogg = 8.0), who 
however included H/3, and Kidder (1991; T^s^ 42500K, 
\ogg — 7.6), who fitted Lyman-a and the Balmer lines H/3 
and H7. 

REJ1036+460, REJ1043+490, REJ1126+183: 

Three more white dwarfs of our observed sample are known 
binaries (RE J1036-f460, RE Jl 126-1-183: Green, Schmidt, 
& Liebert 1986, RE J1043-f-490: Schwartz et al. 1995) and 
show red excesses in our spectra: Since the H/3 lines of the 
white dwarfs are contaminated they were excluded from 
the fit. H7 and higher Balmer lines are virtually uncon- 
taminated. A discussion of these stars and newly discov- 
ered binaries will be given in a forthcoming paper. 

RE J0134-160 (GD 984): Ahhough our spectrum ends 
at only 5600 A a red excess is obvious. Subtraction of 
the theoretical white dwarf flux leaves an M star spec- 
trum. The M dwarf contribution is much smaller than 



in RE J1629-I-780 and we excluded only the cores of H/3 
and H7, filled in by the Balmer line emission, from the fit. 
Bues & Asian (1995) suspected a hot third component in 
REJ0134-160. However, outside of the H/3 and H7 cores, 
the Balmer lines are well reproduced by our best fit with- 
out any indication of a third component. Indeed, a sub- 
dwarf component as suggested by Bues & Asian (1995) is 
almost certainly not present, since its flux would dominate 
in the blue. 

RE J1711+664: Since a late type star is separated from 
the white dwarf by only 2'.'5, we took care to get an uncon- 
taminated white dwarf spectrum. We obtained a spectrum 
under good seeing conditions. No excess is present up to 
the red hmit at 8500 A. 

RE J2207+252: This white dwarf has a red companion 
8'.'5 away. Schwartz et al. (1995) estimated spectral type 
K4V and distance 65 pc from its colors. With the param- 
eters from Table ^ a white dwarf distance of 62 pc results. 
Thus it is likely that both stars form a physical pair. The 
angular separation corresponds to «500AU. 

5.2. Magnetic white dwarfs 

Two white dwarfs of our sample show Zeeman split- 
ting of the Balmer line cores, indicative of a mag- 
netic field. REJ1659-)-440 (PG 1658-^441) is a well 
known star, already analyzed by Schmidt et al. (1992). 
The discovery of the magnetic nature of RE J1440-I-750 
(HS 1440+7518) was announced by Dreizler et al. (1994; 
note the naming confusion corrected in Dreizler et al. 
1995: HS 1412-^6115 should have been HS 1440+7518). 
Although REJ1440+750 was analyzed by V97, they did 
not remark on its magnetic nature. This is likely due to 
the lack of coverage of the Ha line, which displays the most 
pronounced Zeeman effect. 

Flux calibrated spectra of the magnetic white dwarfs 
are displayed in Fig. |. The PC 1658+441 analysis of 
Schmidt et al. (1992) resulted in Teff = 30500 K and log 5 = 
9.35. The magnetic splitting was best reproduced by a 
3.5megagauss (MG) dipole inclined 60° to the line of sight 
(producing a mean surface field strength = 2.3 MG). 
From the linear Zeeman effect we estimated a mean mag- 
netic strength of 8MG for RE J1440+750, consistent with 
the estimate given in Dreizler et al. (1994). 

The temperature and gravity of RE J1440+750 were de- 
rived from a fit of the higher Balmer lines H7 to H77, which 
are less affected by the magnetic splitting. Results are 
given in Table || supplemented by the Schmidt et al. (1992) 
fit of PG 1658+441. Our own fit gave similar results albeit 
with lower accuracy. Our fits of RE J1440+750 can only 
provide a rough estimate of the stellar parameters. Accu- 
rate results can only be expected from a detailed treatment 
of the magnetic effects. 

5.3. Internal, external, and systematic errors 

We have now presented the results of a homogeneous 
analysis of a sample of 46 hot, EUV selected white dwarfs 
based on Balmer line fitting. Three other large samples 
analyzed with the same method were recently published 
by M97, V97, and FKB. Since considerable overlap exists 



Table 2 

RESULTS OF REPEATED OBSERVATIONS 



Observation 1 Observation 2 



RE 


Tcs 


log.g 


Ten 


logg 


1235+233 


46308+837 


7.85+0.08 


46737+670 


7.82+0.06 


1257+220 


39349+317 


7.76+0.04 


38820+159 


7.78+0.02 


1643+411 


28813+120 


8.23+0.02 


28818+111 


8.21+0.02 


1650+403 


37798+245 


7.94+0.04 


39126+413 


8.01+0.05 


1726+583 


52712+927 


8.27+0.06 


54003+669 


8.20+0.05 


2116+735 


50131+384 


7.71+0.03 


54604+906 


7.76+0.06 



NOTES: The error estimates were taken from the fitting 
procedure. 



Table 3 

PARAMETERS OF the MAGNETIC WHITE DWARFS 



RE 


other names 




logs' 


M/Mq 


Bs 


Observation 












Date 


1440+750 


HS 1440+7518 


36154+275 


8.87+0.05 


1.128 


7.7 


10 May 1996 


1659+440^ 


PG 1658+441 


30510+200 


9.36+0.07 


1.311 


2.3 


10 May 1996 



"^Parameters from Schmidt et al. 1992 
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between all four samples, this allows a direct check for sys- 
tematic errors and the individual scatter on a star by star 
basis for white dwarfs hotter than 25000 K. 

Since in contrast to previous works, our analysis is based 
on NLTE model atmospheres, we applied the correction 
vectors given in Fig. || to correct for the LTE assumption. 
Since for the hottest white dwarfs these corrections be- 
come large, while the accuracy of temperature and gravity 
estimates decreases for both LTE and NLTE analyses, one 
should exclude comparison of stars with Teff> 70000 K. 
This does not affect our sample, which has a maximum 
temperature closer to 54000K. 

Differences (after correction to NLTE) in Tofrand logg 
between studies for stars in common with M97, V97, and 
FKB are displayed in Fig. || as function of Teff. The mag- 
netic white dwarfs and the binaries, which show signifi- 
cant contamination of the white dwarf spectrum by the 
companion, are excluded. One can now focus on system- 
atic differences between pairs of studies, e.g. our results 
(NGS) versus FKB (FKB-NGS), V97-NGS, M97-NGS, 
FKB-M97, M97-V97, and so on. However, we chose an- 
other approach and performed an optimization that si- 
multaneously took into account all values from all stars in 
common between any samples. In other words, the val- 
ues given for the systematic differences between our study 
and the M97, V97, and FKB samples form a system for 
direct transformation between, say, FKB and the three 
other samples. The running averages computed this way 
are plotted in Fig. ^ as a solid line. The actual value was 
computed in 1000 K steps for all white dwarfs in the re- 
spective temperature intervals roff±10%. Although this 
curve does not represent the best fit to the data plotted 
in Fig. H alone, it is a fairly good representation of the 
differences computed directly between our measurements 
and those of M97, V97, and FKB. 

Since the distribution is highly non-Gaussian with many 
outliers, as evident in Fig. ^, we decided to adopt an un- 
derlying Lorentzian (or Cauchy) distribution for the opti- 
mization. The tails of the Cauchy distribution are much 
larger than that of the corresponding Gaussian, yielding 
a much lower weight for deviant points (see discussion in 
Press et al. 1986). The dotted lines represent the la confi- 
dence interval of the mean, computed conservatively from 
the rms deviations. 

First, we notice a considerable scatter, larger than ex- 
pected from the internal error estimates (see the discussion 
below). If one ignores the hot end, the agreement between 
the FKB and our temperature scale is good; differences 
are below 1%, smaller than the maximum model differ- 
ences to the Koester models (cf. Sect. 4), which were used 
by FKB. The same atmospheres are used in M97, and it is 
therefore surprising that significant differences with M97 
are present. These trends are most likely caused by dif- 
ferent reduction and analysis techniques. Offsets of the 
same order are found in our comparison with V97, where 
a different LTE model atmosphere code is used. Although 
basically the same input physics is included, this might at 
least partly explain those shifts in Toffa-nd logg. 

Our results are quantified in Tables | and ^ We divided 
the white dwarfs into three groups according to their tem- 
perature: a cool group with Toff < 30000 K, a hotter one 
with 30000 K < Toff < 45000 K, and the hottest consid- 



ered group with 45000 K < Toff < 70000 K. NLTE effects 
are still negligible in the range of effective temperatures 
of the two coolest groups. Mean shifts and the confidence 
range of the mean (computed as described above) are pro- 
vided for these groups and the complete sample (except 
stars with Toff > 70000 K). The shifts discussed above and 
shown in Fig ^ are statistically significant and reach val- 
ues of ~5% in Toff in the hottest bin. In \ogg the differ- 
ence between our results and M97 reach wO.l dex for the 
coolest bin. However, we emphasize that all four analyses 
are based on state-of-the-art model atmospheres and 
fitting techniques. Since there are no strong arguments 
to favor or discard one analysis, it seems these are the 
systematic shifts characteristic of modern analyses of hot 
white dwarfs. 

If we take the systematic shifts into account, we can 
use the samples to derive estimates of the observational 
scatter, which can be compared with the internal error es- 
timates. Since we have approximately the same scatter for 
different combinations of the samples, we compute a mean 
scatter cr^ig for all possible combinations (weighted by the 
number of stars in common) . For this purpose we correct 
for the systematic shifts calculated for each of the three 
Toff bins, whose temperature intervals are given in Table | 
The individual measurement errors Cind add quadratically, 
and if we assume that the inherent scatter is the same for 
all analyses (and we found no reason to discard this as- 
sumption), the individual measurement errors can be es- 
timated as CTind — -^fdiff. Not surprisingly, the smallest 
scatter is found for the "cool" group (TefT < 30000 K) with 
o'(Tofr) — 2.3% and cr(logg) — 0.07 dex. It increases to 
o'(7off) = 3.3% and cr(logg) = 0.13 dex for the hottest 
bin. This trend is expected from the statistical analysis 
presented in FKB (their Fig. 1). However, the values are 
larger by a factor of three or more than the internal param- 
eter errors for a well exposed spectrum (see e.g. Table p. 
Therefore, we conclude that the accuracy is not limited 
by the noise for good spectra, and we suggest that other 
effects, such as details of the extraction or fluxing and 
normalization procedures, contribute more. Considering 
these systematic uncertainties, the 0.3 dex difference be- 
tween the gravity determinations of FKB and Napiwotzki 
et al. (1993) for HZ 43 A is only a 1.5a deviation and there- 
fore not as serious as considered by FKB. 



5.4. Mass distribution 
5.4.1. Derivation and Sample Comparisons 

Once the temperature and gravity of the white dwarfs 
are known the mass can be determined from theoretical 
mass-radius relations. The recent investigations of M97, 
V97, and FKB based their interpretation on the model 
sequences of Wood (1995). The applied models have H- 
and He-layer masses of lO^^AfyvD and IO^^A/wd, respec- 
tively ("thick layers") and a carbon core. These model 
sequences are based on a set of pre-white dwarf models 
with end masses of O.GMq, O.SMq and 0.95Mq computed 
by Kawaler (see Wood 1994). Starting models for other 
masses were constructed by homology transformations. 

The mass-radius relations we use are based on the evo- 
lutionary calculations of Blocker (1995). These models 
are calculated ab initio from the main sequence. The 



9 



Table 4 

SYSTEMATIC TEMPERATURE DIFFERENCES 





M97-NGS 

A Teff /Teff 


V97-NGS 

A Teff /Teff 


FKB-NGS 

A Teff /Teff 


Teff< 30000 K 
30000 < Teff < 45000 K 
45000 < Teff < 70000 K 
all with Teff < 70000 K 


-0.023±0.006 
-0.013±0.003 
0.036±0.006 
-0.006±0.003 


0.032±0.005 
0.017±0.003 
0.050±0.006 
0.024±0.003 


0.007±0.005 
0.004±0.003 
0.046±0.007 
0.012±0.003 



NOTES: For each temperature range, each column gives the average 
difference over the respective interval, and the confidence interval of the 
average difference. 



Table 5 

SYSTEMATIC GRAVITY DIFFERENCES 







M97-NGS 
Alog.g 


V97-NGS 
A log 5 


FKB-NGS 
Alogg 


Teff < 30000 K 
30000 < Teff < 45000 K 
45000 < Teff < 70000 K 
all with Teff < 70000 K 


-0.112±0.017 
-0.086±0.017 
~0.066±0.030 
-0.084±0.013 


-0.001±0.016 
0.026±0.016 
-0.067±0.026 
-0.002±0.012 


-0.064±0.018 
-0.034±0.016 
-0.062±0.027 
-0.047±0.012 



NOTES: For details see Table | 



Table 6 

MEAN SCATTER a-.^A of TEMPERATURE and GRAVITY 
DETERMINATION 



Teff < 30000 K 0.023 0.075 

30000 < Teff < 45000 K 0.023 0.101 

45000 < Teff < 70000 K 0.034 0.128 

all with Teff < 70000K 0.026 0.107 
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resulting white dwarfs show important differences com- 
pared to Wood's models. In contrast with the carbon 
core models widely used for mass determination, the de- 
generate core contains a mass-dependent mixture of car- 
bon and oxygen (and trace elements). However, note that 
Wood (1995) also presented some white dwarf sequences 
with C/0 cores. The H- and He-layer masses resulting 
from evolutionary calculations depend on the stellar mass 
(Blocker 1995; Blocker et al. 1997). While the canonical 
thick hydrogen layer mass of 1O~'*M0 is met for a O.GM© 
white dwarf, it decreases from several 10~^Mq for O.SM© 
down to 1O~^M0 for l.OM©. This results in a mass-radius 
relation which is steeper than a relation based on a con- 
stant layer mass. Consequently, the derived mass distri- 
bution will be narrower, if evolutionary layer masses are 
used. 

The lower mass limit for a C/0 core white dwarf is 
0.46Mq, the limiting mass for central helium burning in 
low mass stars (Sweigart et al. 1990). Thus, white dwarfs 
with a lower mass possess a helium core. At the current age 
of the universe He white dwarfs have not been produced 
by single star evolution, but are the result of binary evolu- 
tion where the hydrogen-rich envelope was stripped away 
along the (first) red giant branch (Kippenhahn 1967; Iben 
& Tutukov 1986). We use the recent evolutionary models 
of He white dwarfs calculated by Driebe et al. (1998) based 
on this scenario. 

The combined tracks of Blocker (1995) and Driebe et 
al. (1998) cover the mass range O.I8M0 to 0.94Mq. We 
supplemented this set with the 1.0, 1.1, and 1.2Mq car- 
bon core sequences of Wood (1995) with "thin" layers. 
Since the hydrogen envelope mass decreases with increas- 
ing white dwarf mass in Blocker's models, their highest 
mass models effectively correspond to "thin layer" mod- 
els. Since high mass white dwarfs are already close to the 
zero-temperature configuration, the remaining departures 
from consistent evolutionary calculations can be ignored 
for practical purposes. 

The position of the analyzed white dwarfs in the temper- 
ature/gravity plane is shown in Fig. ^ along with the tracks 
used for the mass determination. The individual white 
dwarf masses are given in Tab le pi , and the resulting mass 
distribution is shown in Fig. OT We have also redeter- 
mined masses for the BSL sample with the Blocker/Driebe 
mass-radius relations and supplement Fig. |l^ with the re- 
sulting distribution. 

Our mass distribution possesses the same basic features 
as found in the EUV selected samples of M97, V97, and 
FKB. The sharp peak centered at w O.59M0 is in principal 
agreement with the earlier investigations of KSW, Weide- 
mann & Koester (1984), and BSL, with a sharp falloff 
towards lower masses, and a less steep decline towards 
higher masses with a tail of high mass white dwarfs. Con- 
sidering the relatively large observational scatter we found 
in Sect. |5.3| , the underlying distribution may be extremely 
sharp. The mean mass is O.67M0. However, this value is 
strongly biased by the few very high mass white dwarfs. 

We decided to follow the recipe of FKB and fitted the 
mass peak with a Gaussian. Although a Gaussian is not 
a good representation of the white dwarf mass distribu- 
tion, it gives a robust estimate of the peak mass. The fit 
result is, in principle, dependent on the adopted binning 



of the mass intervals. However, in some test calculations 
we found effects exceeding a few O.OOIM© only if the bin 
width was larger than O.O5M0. The values given in this 
paper were derived by fitting Gaussians to virtually un- 
binned mass distributions (formal bin width O.OOIM©). 

We were concerned by the possibility that the different 
fraction of low and high mass white dwarfs found in differ- 
ent samples (see discussion below) may skew the Gaussian 
fit of the main peak to higher or lower masses. We tested 
this by fitting the BSL, and FKB and our mass distribu- 
tions with multi-component Gaussians, which fitted the 
secondary high and low mass peaks separately. We found 
no deviation higher than O.OOSA/© in any case and con- 
cluded that the single Gaussian fit of the main peak yields 
a rather stable estimate. 

With this method, we derived a peak mass of our sam- 
ple of O.589M0. We reanalysed the FKB, M97, and V97 
samples with our mass-radius relations. We applied the 
corrections to NLTE and excluded the white dwarfs with 
temperatures in excess of 70000 K. The peaks masses are 
O.555M0 (FKB), O.535Af0 (M97), and 0.582 (V97). These 
differences reflect the systematic differences in Toff and 
log 5 determinations discussed above and probably differ- 
ent selection criteria in the case of FKB. 

The probable He core white dwarf RE JOS 12-004 is the 
only object in our sample with a mass below 0.5Mq, but 
four white dwarfs have masses in excess of I.OMq. The fre- 
quency of low and high mass white dwarfs found in our and 
other EUV selected samples is qualitatively different from 
that found in optically selected samples such as BSL. BSL 
detected in their sample of 129 white dwarfs 16 objects 
with masses below the limiting mass for C/0 core white 
dwarfs and a total of 28 white dwarfs with M < 0.5Mq. 
These numbers are reduced to 10 and 16 respectively, if 
we redetermine the masses with the mass-radius relations 
of Blocker (1995) and Driebe et al. (1998). However, that 
is still a much higher fraction than we find in our EUV 
selected sample. 

5.4.2. Cooling Rates and Detection Probability 

Why do optical and EUV selected samples contain dif- 
ferent fractions of high and low mass white dwarfs? All 
helium white dwarfs from the BSL sample have tempera- 
tures below 30000 K; the only helium white dwarf candi- 
date in our sample has Tqb «32000K. This implies that 
the detection probability of low mass white dwarfs in EUV 
selected samples is much lower, because the fraction of low 
mass white dwarfs at higher T^e is much lower than it is for 
Tcff < 30000 K. This viewpoint is supported by the helium 
white dwarf sequences calculated by Driebe et al. (1998). 
We have plotted the variation of temperatures with time 
for white dwarfs with helium and C/0 cores in Fig. ^ 
The detection probability in a given evolutionary stage 
goes as the inverse of the rate of temperature decrease at 
that stage. New He white dwarfs cool down rapidly to 
a certain (mass dependent) temperature, whereafter the 
cooling rate drops drastically and the star remains on a 
temperature plateau for a long time. This behavior is 
most pronounced for the lowest masses. Hydrogen shell 
burning plays an important role for this plateau phase. Al- 
though the energy production by hydrogen burning drops 
down dramatically after the star enters the white dwarf 
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sequence, it still produces a significant fraction of the to- 
tal luminosity, enough to bring the cooling to a standstill. 
We can conclude from Fig. ^ that it is extremely unlikely 
to detect a helium white dwarf during its first rapid cool- 
ing phase. It is easier to find them during their plateau 
phase or afterwards. This means there is a strong bias 
towards detecting helium white dwarfs in samples of rel- 
atively "cool" DA stars (like BSL) instead of the overall 
much hotter EUV selected samples. 

We can explain in a similar fashion why high mass white 
dwarfs are found preferentially in the 30000 K < T^s < 
50000 K range. Due to the drop in neutrino cooling, the 
evolutionary speeds of massive white dwarfs are relatively 
low in this range, as demonstrated by the O.94OM0 track, 
the most massive stellar remnant calculated by Blocker 
(1995). This increases the detection probability of more 
massive white dwarfs in EUV selected samples. 



5.4.3. Interstellar Absorption 

Even if we restrict our comparisons to hot subsamples of 
white dwarfs, the fraction of massive white dwarfs is much 
higher in the EUV selected samples than it is in optically 
selected ones. An extensive discussion of this problem is 
carried out by FKB. FKB argue that the number of high 
mass white dwarfs discovered per sky area is similar for 
the whole sky EUV surveys and the optical Palomar-Green 
(PG) survey (Green et al. 1986), which covered roughly 
25% of the sky. In other words, a whole sky version of 
the PG survey should have discovered all but one of the 
EUV detected white dwarfs with M > I.IMq. At a given 
temperature, white dwarfs with lower masses are larger, 
consequently more luminous, and detectable to larger dis- 
tances. For instance at 30000K, a white dwarf of O.5M0 is 
detectable at distances ^ 30% larger than a 0.7Mq white 
dwarf. At the maximum distance sampled by, e.g., the PG 
survey, the probability is very high that interstellar matter 
effectively absorbs all EUV radiation. Thus the dominant 
effect in EUV surveys is not a selection for massive white 
dwarfs as suggested by V97, but a selection against low 
mass white dwarfs, due to a sample volume strongly af- 
fected by interstellar absorption. 



5.4.4. Temperature Dependence of the Derived Mass Peak 

FKB noted a moderate apparent trend of the white 
dwarf peak mass with temperature. Since their sample 
also included optically selected white dwarfs, they ana- 
lyzed a considerably higher fraction of stars with Tcff < 
25000 K and therefore have a larger temperature base- 
line than pure EUV samples. Analogous to this find- 
ing, V97 found an O.O3M0 offset between the peak of 
the mass distribution in their EUV selected sample and 
the (cooler) white dwarfs analyzed by BSL. In both cases 
this temperature dependence decreases if Wood's mod- 
els with "very thin" layers (Mhb = 10~'''Mwd, no hy- 
drogen layer) instead of the canonical thick layer models 
(Mne = 10-2Mwd,Mh = IO'^Mwd) are used. How- 
ever, part of the trend noted in V97 may be related to 
the systematic differences between the analyses discussed 
in Sect 5.3. Since V97 compared results from two dif- 
ferent analyses, it is not possible to judge if the offset is 



real or just an artifact. The case of FKB is stronger, be- 
cause they based their case on a homogeneously analyzed 
sample. FKB warned that the temperature dependence of 
the mass distribution peak may be caused by inadequa- 
cies remaining in the model atmospheres. A more detailed 
discussion is given in the next section. 

Some evidence that part of the temperature trend stems 
from the analyses comes from our intercomparisons pre- 
sented in Sect 5.3. We showed that systematic differ- 
ences between the four investigated samples exist that vary 
with effective temperature. This might mimic a temper- 
ature dependence of the sample peak mass. A gravity 
offset of 0.1 dex transforms into mass offsets of O.O5OM0, 
0.044Mq, and O.O34M0 for a 0.6Mq white dwarf with 
25000 K, 40000 K, and 60000 K, respectively. Even a rela- 
tively small systematic \ogg difference of 0.05 dex (cf. Ta- 
ble ||) corresponds to a O.O2M0 offset. The scatter a of 
individual gravity determinations reported in Table ^ cor- 
responds to a{M) w O.O4M0, nearly independent of Teg. 



5.4.5. Cooling Tracks 

Another part of the T^s dependence of the mass peak 
may be caused by the use of Wood's cooling tracks by FKB 
and V97. These tracks are not completely self-consistent 
since not all model sequences were not calculated ab ini- 
tio from the main sequence. This is not important for 
cooler white dwarfs but may cause deviations from ab ini- 
tio tracks especially for the hottest white dwarfs, where 
structure depends sensitively on the evolutionary history 
(Blocker & Schonberner 1990). 

We redetermined the masses of the BSL white dwarfs 
with the Blocker/Driebe mass-radius relations and derived 
a peak mass of O.559M0 . That is O.O3M0 lower than our 
peak mass, but having in mind the range of mass determi- 
nation derived for the EUVE selected samples, we cannot 
consider this a significant difference. FKB divided their 
sample into cool {T^s < 35000 K) and a hot (35000 K 
< Toff <75000K) subsamples and derived a O.O29M0 
higher peak mass for the hot sample. Our reanalysis of 
the FKB LTE results with the Blocker/Driebe mass-radius 
relation yields virtually the same offset: O.O25Af0. How- 
ever, the difference is brought down to O.OIOM0 (O.55OM0 
vs. O.56OM0) if we apply corrections to NLTE. One can 
imagine that a difference of this order (if significant at all) 
can easily be produced by our neglect of metallicity effects 
(cf. Lanz et al. 1996, Barstow et al. 1998). 

Therefore, given the combination of sample selection 
and systematic effects in analyses to date, our results do 
not confirm the presence of intrinsic, systematic mass dif- 
ferences between hot and cool white dwarfs. 



6. CONCLUSIONS 

We have obtained temperatures, gravities, and masses 
for a sample of 46 extreme ultraviolet selected DA white 
dwarfs. These data complement a near infrared survey for 
low mass companions. The stellar parameters were deter- 
mined by fitting the hydrogen Balmer line profiles with 
NLTE model spectra. A map of LTE correction vectors 
was constructed, which allows the transformation of LTE 
results to the NLTE scale. These shifts become impor- 
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tant for temperatures above 50000 K (see also Napiwotzki 
1997). 

Three recent analyses of DA white dwarfs have a consid- 
erable overlap with our sample and allow a direct check for 
systematic errors and individual measurement scatter. Al- 
though all four analyses applied similar model atmospheres 
and fitting techniques, we recognized systematic shifts up 
to 5% in temperature and 0.1 dex in gravity determina- 
tions. These systematic errors have to be considered when 
the results arc interpreted and compared to other studies. 
The individual measurement errors increase with tempera- 
ture: cr(Teff) from 2.3% to 3.3% and a{\ogg) from 0.07 dex 
to 0.13 dex. This is larger than the typical errors computed 
with fit procedures for a well-exposed spectrum and in- 
dicates that accuracy is limited by noise in combination 
with other effects such as details of the extraction, fluxing 
or normalization procedures. It is unclear if these a val- 
ues indicate lower limits, or scale (at least partially) with 
the S/N ratio of the spectra. It will be a challenge to re- 
solve this question. Repeated observations performed by 
one observer, at one telescope, reduced the same way and 
analysed with one method as performed by BSL, FKB, 
and us (Table 2) do not reveal all effects. Independent 
observations analysed independently are necessary. 

Masses have been inferred from theoretical mass-radius 
relations based on the evolutionary calculations of Blocker 
(1995) for C/0 white dwarfs and Driebe et al. (1998) for 
He white dwarfs. An important feature of these calcu- 
lations are the hydrogen and helium layer masses which 
depend on stellar mass. We find a sharp peak centered at 
SSO.59M0 in agreement with the previous investigations 
of KSW and BSL. We redetermined masses of the BSL 
white dwarfs with the Blocker/Driebe mass-radius rela- 
tion and derived a peak mass of O.56M0. At first glance 
this seems to confirm the systematic offset reported by 
V97. However, this offset could be explained by system- 



atic differences of the model atmosphere analyses as well. 
Our reanalysis of the homogeneous FKB sample showed 
that the temperature dependence of the mass peak nearly 
vanishes, when corrected for their LTE assumption. We 
conclude that the observational data presented here and 
in similar studies are well explained by canonical stellar 
evolution theory, i.e. white dwarfs with thick envelopes. 

We find only one object (RE J0512-004) with a mass be- 
low O.5Af0, which is a possible helium core white dwarf, 
but four white dwarfs with masses in excess of I.OMq. 
This ratio of high- and low-mass white dwarfs is quite dif- 
ferent from that found by BSL. This can partly be under- 
stood as the result of variable evolutionary time scales of 
high and low mass white dwarfs. We agree with FKB that 
another part can be explained by a selection against low 
mass white dwarfs in EUV selected samples. 

The principal aim of our project is the search for and 
study of cool main sequence companions of our white dwarf 
sample. Five white dwarfs already show red contamina- 
tion in our optical spectra. A (•oniprehcnsive discussion 
of these and other binaries will be given in a forthcoming 
paper presenting IR photometry for the white dwarf sam- 
ple. In this article we provided the basic white dwarf data 
necessary to interpret our binary sample. 

We thank T. Driebe, F. Herwig, and T. Blocker for pro- 
viding us with their tracks and computing the white dwarf 
masses and D. Kocster, who made some model spectra 
available for our model comparison. We are grateful to 
Steward Observatory for an unexpected, generous award 
of six nights of unclaimed 2.3-m time, which made possible 
simultaneous optical and infrared observations of our sam- 
ple of stars. We thank Perry Berlind for obtaining several 
white dwarf spectra. PJG acknowledges support through 
NASA Contract NAS8-39073 (ASC). 
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Fig. 1. — Spectra of the white dwarfs observed at Steward observatory 
Fig. 2.— Spectra of the MMT sample 

Fig. 3. — LTE offsets. The differences are magnified three times. The vectors give the correction which must be apphed 
to transform LTE results to the NLTE scale. 

Fig. 4. — Balmer line fits for a representative set of white dwarfs. 

Fig. 5.— Flux calibrated spectrum of the DA+M pairs RE J0134-160, RE 1036+460, RE J1043+490, RE J1126+183, and 
RE J1629+780. From top to bottom we show the composite spectrum, the best fit of the white dwarf, and the resulting 
spectrum of the M star. The spectrum of the M star companions of RE J0134-160 and RE J1126+183 are multiplied by 
3. 

Fig. 6. — Spectra of the magnetic white dwarfs PG 1658+441 and REJ1440+750. The position of the a components of 
the Balmer lines are indicated. 

Fig. 7. — Best fit of the magnetic white dwarf RE J1440+750. As indicated by the dashed line, H/3 was not used for the 

parameter determination. 

Fig. 8. — Differences in temperature (left), and gravity (right) between the M97, V97, and FKB samples and our own on 
a star by star basis. The same comparison is carried out for the LTE studies of M97 and V97 in the bottom panel. The 

smoothed average of the differences is plotted with solid lines, while the standard error of the mean difference is indicated 
by the dotted lines. The T^s values used for the x-axis are NLTE from this paper (NGS), except for the bottom panel, 
which uses V97 (corrected to NLTE). 

Fig. 9. — Effective temperature and gravity of our white dwarf sample compared with evolutionary tracks. Solid lines are 
Blocker (1995); dashed lines are Driebe et al. (1998); dotted lines Wood (1995). The magnetic white dwarfs are marked 

by open symbols. 

Fig. 10. — Mass distribution of white dwarfs (binned over 0.05 Mq). The solid lined histogram indicates the result for 
our sample. The dotted line is the BSL sample with masses redetermined from the Blocker/Driebe relations (normalized 
to the same sample size). The Gaussian curve represents the best fit to our mass distribution as described in the text. 

Fig. 11. — Time evolution of stellar surface temperature for helium core white dwarfs (tracks from Driebe ct al. 1998; solid 
lines) and C/0 white dwarfs (tracks from Blocker 1995; dashed lines). The tracks are labeled with the stellar masses. 
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